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Ground State Entropy in Potts Antiferromagnets and Chromatic 
Polynomials 

R. Shrock a * and S.-H. Tsai a * 
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We discuss recent results on ground state entropy in Potts antiferromagnets and connections with chromatic 
polynomials. These include rigorous lower and upper bounds, Monte Carlo measurements, large-q series, exact 
solutions, and studies of analytic properties. Some related results on Fisher zeros of Potts models are also 
mentioned. 
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1. Introduction 

Nonzero ground state entropy, So ^ 0, is an 
important subject in statistical mechanics, as an 
exception to the third law of thermodynamics 
(e.g. A physical example is ice. A simple 

model exhibiting ground state entropy is the q- 
state Potts antiferromagnet (PAF) on a lattice, 
or more generally a graph G, for sufficiently large 
q. The zero-temperature partition function satis- 
fies 



Z{G,q,T = 0) PAF = P(G,q) 



(1) 



where P(G, q) is the chromatic polynomial ex- 
pressing the number of ways of coloring the ver- 
tices of G with q colors such that no two adja- 
cent vertices have the same color. Let {G} = 
linin^oo G. The ground state degeneracy per site 
W, given by So = fcs In IT, is 



W({G},q) = km P(G,q) 



l/n 



(2) 



There are very few nontrivial exact solutions for 
W . We have obtained a number of new results on 
W , including rigorous bounds, large-g series ex- 
pansions, Monte Carlo measurements, exact solu- 
tions, and studies of analytic properties P|-pO[. 
Besides physical (positive integral) q, it is of in- 
terest to consider IT as a function of complex q; 
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W is an analytic function of q except on a contin- 
uous locus B (B may be null, and there may also 
be isolated singularities of W). As n — > oo, the 
locus B forms via coalescence of a subset of zeros 
of P(G, q). W is determined via ^ in the region 
R\ reached by analytic continuation from the real 
q axis for q > x(G), where x(G) is the chromatic 
number of G, i.e., the minimum number of colors 
needed to color G with the above constraint. In 
other regions separated from R\ by nonanalytic 
boundaries comprising B, only \W\ can be deter- 
mined. There is a subtlety in the definition of W, 
since at certain special points one encounters the 
noncommutativity of limits [3j 

I 1 /" ^ Km Km P(G,g) 1/ 



lim lim P(G,qY 

n — >oo q — >q s 



^ lim 



(3) 



At such points, we use the second order of limits 
to define W. Physically, one finds W > 1, i.e., 
So > for q > x(G). 

2. Bounds, Series, and Monte Carlo Mea- 
surements 

We have proved rigorous lower and upper 
bounds on W for a number of lattices A 
As an example, for the honeycomb (he) lattice, 
we get 

W{hc, q) t < W(hc, q) < W(hc, q) u (4) 
where 

(^-Sg^lO^-lOg + S) 1 / 2 
W(hc, q)i = — (5) 
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and 

W(hc,q) u = (q 2 -3q + 3) 1/2 (6) 

These bounds are very restrictive even for moder- 
ate q, as is clear from the fact that the first three 
terms in a large-g expansion coincide. Although 
a bound on a given function need not, a priori, 
coincide with a series expansion of that function, 
we find that the lower bound (^) coincides with 
the first eleven terms of the large-g expansion for 
W(hc,q). 

Accordingly, we have extended our study of 
the lower bound and have discovered and proved 
a generalization applicable to the full set of 
Archimedean lattices ||. An Archimedean lat- 
tice is a uniform tiling of the plane with one or 
more regular polygons such that all vertices are 
equivalent to each other. It can be specified by 
the ordered sequence of polygons Pi traversed by 
a circuit around any vertex: A = fl^p"*. -Let 
^2 a,i = a,i tS and Vi = ai^jpi. Then our general 
lower bound is 

w((U P ru)>^^ (7) 

i 

with 

D k (q)=J2(-iy ( k ~% k - 2 - s (8) 

s=0 ^ ' 

We have calculated large-g series expansions for 
a number of Archimedean lattices and have com- 
pared the lower bounds with these series For 
the square, triangular, and honeycomb lattices we 
have carried out Monte Carlo measurements of 
W(A, q) for q values up to 10 and have found 
that even for moderate values of q, the upper and 
lower bounds bracket the measured values quite 
closely 0. 

3. Analytic Structure of W({G},q) 

We have calculated exact solutions for 
Vy({G}, q) for a number of families of graphs and 
have studied their analytic structure. A general 
form for P(G, q) is 

P(G,q)=c (q) + J2 c M( a M) tjn (9) 



where Cj(q) and aj(q) are certain functions of q. 
Here the a,j(q) and Cj^o(q) are independent of n, 
while Co (q) may contain n-dependent terms, such 
as (—1)™, but does not grow with n like (const.) 71 
with | cons£.| > 1. A term ai(q) is leading if it 
dominates the n — > oo limit of P(G, q); in partic- 
ular, if N a > 2, then it satisfies |a£(g)| > 1 and 
\a t (q)\ > \aj(q)\ for j ^ i, so that \W\ = \a e \^ . 
The locus B occurs where there is a nonanalytic 
change in W as the leading terms ai in eq. (^) 
changes. For a given {G} one can ask various 
questions about B: (i) does it separate the q 
plane into separate regions? (ii) is it compact 
or noncompact? (iii) how many disjoint compo- 
nents does it contain? (iv) does it have multiple 
points where several branches cross? (v) does it 
cross the positive real axis or contain a segment 
lying along this axis, and if so, what is the maxi- 
mal real value, q c , in £>? We have answered these 
questions for various families P,p|.|(|[7||To"[ . 

For example, such families as circuit and lad- 
der graphs yield B that do separate the q plane 
into various regions. With M. Rocek, we have 
used a generating functional method to obtain 
exact solutions for W on infinitely long, finite- 
width homogeneous strip graphs ||. For these 
we find that the loci B consist of arcs that do not 
separate regions. 

In general, the zeros of P(G, q) only merge to 
form B in the n — > oo limit, but for an inter- 
esting set of families, they actually lie exactly 
on B even for finite n. We have constructed 
an infinite set of these families, called p-wheels, 
(Wh)n^ = K p + C n -p, where K r is the "com- 
plete" graph (each of whose vertices is connected 
to all others), C r is the circuit graph, and G + H 
connotes joining of vertices of G to those of H . 
We proved that the zeros ofP((Wh) { P\q) lie pre- 
cisely on B, which is the circle \q— (p+l)\ = 10. 
This is reminiscent of the Yang-Lee circle theorem 
|nj although different in certain respects. 

For {G} with noncompact B passing through 
z = 1/q = 0, no large-g series expansions exist. 
Since these series are very useful for regular lat- 
tices, it is important to understand the properties 
of a family {G} that yield a noncompact B. We 
have given a general condition for B to be non- 



3 



compact and have constructed a number of these 
families 0; one simple infinite set of families is 
G = (Kp)b + H, where b signifies the removal of 
b bonds from a vertex of K p . For example, for 
1 < b < p — 1 and H = T r (the tree graph on r 
vertices), in the limit r — > oo, B is the vertical line 
with Re(q) = p+^, or equivalently, in the z = 1/q 
plane, the circle \z — z c \ — z c with z c = 2/(2p+l). 
Additional families can be generated by various 
types of homcomorphic expansions of the basic 
G = (K p )b + H (homcomorphic expansion = in- 
sertion of degree-2 vertex on a bond). 

4. Approach to 2D thermodynamic limit 

Using exact solutions of W for infinitely long 
strip graphs of finite width L y , we have ex- 
plored the approach of W to its 2D thermody- 
namic limit ||. We showed that the approach 
of W to its 2D thermodynamic limit as L y in- 
creases is quite rapid; for moderate values of q 
and L y ~ 4, W{A, (L x — oo) x Ly.q) is within 
about 0(lCr 3 ) of the 2D value W(A, (L x = 
oo) x (Ly = oo), q) for periodic transverse bound- 
ary conditions (b.c). The approach of W to the 
2D thermodynamic limit was proved to be mono- 
tonic (non-monotonic) for free (periodic) trans- 
verse b.c. An application to compute central 
charges for cases with critical ground states was 
noted. 

5. Complex-temperature properties 

We have calculated Fisher zeros of the partition 
function for the Potts model on the square lattice 
for several q values and related the inferred CT 
phase boundaries to locations of singularities in 
the thermodynamic functions obtained from anal- 
yses of low-temperature series expansions [fl2f . 
These studies have been extended to honeycomb, 
triangular, kagome, and diced lattices |p|[l4|]. Al- 
though in general CT singularities have rather 
different properties than physical critical singu- 
larities in spin models, such as violation of uni- 
versality [lq] , we have found an exact mapping 
that relates the free energy of the q-state Potts 
antifcrromagnet on a lattice A for the full tem- 
perature interval < T < oo and the free energy 



of the (/-state Potts model on the dual lattice for 
a semi-infinite CT interval, — oo < ad < — (q — 1), 
where a,j = (a+q— l)/(a— 1) and a = e K . Hence, 
for this interval, CT singularities of the free en- 
ergy are equivalent to physical critical singular- 
ities. Effects of next-nearest-neighbor couplings 
were studied in Ref . jl6| . 
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